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To understand newly discovered superconductivity in Fe-based systems, we investigate electronic 
structure and magnetic properties of Fei+a;Te using first-principles density functional calculations. 
While the undoped FeTe has the same Fermi surface nested at (tt, tt) as in Fe arsenides, doping 
by ~ 0.5 electrons reveals a novel square-type Fermi surface showing a strong (7r,0) nesting and 
leading to a different magnetic structure. Our result strongly supports the same mechanism of 
superconductivity in chalcogenides as in the arsenides, reconciling theory to existing experiments. 
Calculated magnetic interactions are found to be critically dependent on doping and notably different 
from the arsenides. 

PACS numbers: 74.70.-b, 71.18.+y, 71.20.-b, 75.25.+Z 



Shortly after the discovery of a novel high tempera- 
ture superconductor, LaFeAsOi_a;F^ with Tc ^27 K 
many different types of iron-based superconductors have 
been reported. Now the highest Tc reaches up to ~55 
K 0, and there are four different structural classes: so- 
called 122-0, and ll-typeQ structures besides 
the original 1111 type. Although tremendous research ac- 
tivities devoted to this field over the year have shed light 
on their intriguing physical properties, our understand- 
ing of superconductivity here and its interplay with mag- 
netism is still far from being complete. One of the most 
important properties which was found in these systems 
is the Fermi surface nesting whose nesting vector corre- 
sponds to the antiferromagnetic (AFM) ordering vector 
of the undoped magnetic phase [1, 0, H, 0] . Density Func- 
tional Theory (DFT) calculations show that all of the 
four classes of these materials share this common feature 
in their electronic structure 0, lid, [ul, which strongly 
suggests the superconductivity is exotic and is mediated 
by spin fluctuations 

Along this line, one of the most interesting questions 
arises in the 11 -type Fe chalcogenide family: Fe(S,Se,Te) 
d, US, El, Q [H [M [13, m • In spite of their same crys- 
tal structure represented by the 2-dimensional Fe square 
lattice and the same Fermi surface nesting predicted by 
DFT calculations fll)|, Fe chalcogenide superconductors 
exhibit notable differences from the arsenides. The mag- 
netic structure found in their parent compound, Fei+^jTe, 
strongly tackles the spin fluctuation theory because the 
experimentally observed magnetic ordering is fairly dif- 
ferent from that of parent arsenide compounds. Although 
DFT calculations predict the same Fermi surface topol- 
ogy, a recent neutron experiment [l9l | shows that, at a 
small X ~ 0.068, Fei+^^Te has a rotated and double stripe 
AFM order. Importantly, this magnetic structure cannot 
be matched with the (tt, tt) nesting as found in all the ar- 
senide materials and predicted by previous calculations, 
but requires (tt, 0) nesting which has never been reported. 
This 'missing nesting' remains as a puzzle in the study 
of Fe-based superconductors. Therefore it is not a big 
surprise that some papers speculate about a different su- 
perconducting mechanism for Fe chalcogenides from the 



arsenides, and the reinvestigation of the electronic struc- 
ture and magnetic properties for Fei+a;Te is of crucial 
importance [lH, |20, 1211, [22| • 

To address these issues, we study electronic struc- 
ture and magnetic interactions of Fei+2,Te using first- 
principles DFT calculations. Since it is noted experi- 
mentally [l^ that FeTe has always some amount of ex- 
cess Fe atoms we perform doping dependent calculations 
to understand their effect on the electronic and magnetic 
properties. Our results show that Fci+^Te has a different 
Fermi surface topology as a function of doping, and even- 
tually a novel (tt, 0) nesting appears at a doping level of 
(5 ~ 0.5 electrons while the (7r,7r) nesting is largely sup- 
pressed. This (tt, 0) nesting exactly matches with the 
double-stripe AFM order found in neutron experiments. 
Our result strongly supports the spin fluctuation medi- 
ated superconductivity for Fe chalcogenides, reconciling 
theory to existing experiments which showed significant 
differences between pnictides and chalcogenides. It is 
also found that magnetic interactions depend on dop- 
ing, and, at the same level of doping 5 ~ 0.5, the cal- 
culated exchange couplings become consistent with the 
double stripe phase. While the second nearest neighbor 
AFM couphng (J2a) is strongest, the first (Jib) and sec- 
ond neighbor ferromagnetic (FM) interactions (J2h) are 
also significant, which is different from arsenides. 

We performed the first-principles DFT calculations 
within local density approximation (LDA) for the 
exchange-correlation energy functionals [11]. The full 
potential linearized-muffin-tin-orbital (LMTO) method 
has been used [1^ . To estimate the exchange interaction 
strengths between Fe moments, we performed a linear 
response calculation [H, [H, 113], which has been suc- 
cessfully applied earlier to Sd transition-metal oxides, 5/ 
metallic alloys [13, HI], and the other Fe arsenides [29| . 
To simulate electron doping a rigid-band approximation 
is utilized. The calculation of band dispersions, Fermi 
surfaces, and Stoner response functions have been done 
with the non-magnetic unit cell (dotted squares in Fig.[T|) 
and using experimental lattice parameters as in the pre- 
vious study ,lli]. To estimate the exchange constants we 
performed a spin-polarized calculation with an enlarged 
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FIG. 1: (Color online) The schematic picture of spin structure 
for (a) iron arsenides parent materials and (b) FeTe. The 
dotted squares correspond to the non-magnetic unit cell and 
the arrows represent the spin directions. The Jia, Jib, Jia, 
and J^b refer to the nearest neighbor AFM, nearest neighbor 
FM, next nearest neighbor AFM, and next nearest neighbor 
FM exchange interaction, respectively. 
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unit cell containing four Fe atoms, and the experimen- 
tal 2:(Te) was used which well reproduces the observed 
moment and is consistent with our previous study of Fe 
arsenides [2^. 

Fig. [T] summarizes spin structures found in the par- 
ent materials of arsenide superconductors (Fig.[T]^a)) and 
Fei+^Te for small x (Fig.fHb)) In Fei+^Te, the FM 
spin stripes are doubled and rotated by 45 degree with 
respect to the non-magnetic unit cell (smallest dotted 
squares). According to the neutron scattering experi- 
ment by Li et al. [iQj, this doubled stripe phase realizes 
at the smallest possible x ~ 0.068, i.e. close to the sto- 
ichiometric FeTe (x — 0). The Jia, Jib, J2a, and J^h 
represent the first nearest neighbor AFM, FM, second 
nearest AFM, and FM interactions, respectively. While 
the second neighbor coupling is always AFM in arsenides 
(Fig. Hlja)), both FM and AFM second neighbor cou- 
plings exist in FeTe (Fig. [Hb)). As x increases further, 
the spin ordering becomes incommensurate at a: ~ 0.141 
and the incommensurate ordering vector depends on x. 
Importantly this different magnetic structure found in 
Fei+xTe cannot be matched with the (tt, tt) nesting which 
is common for AFM parent materials of all the arsenide 
superconductors. The different spin structure found in 
Fei+a;Te remains as a puzzle because DFT calculation 
predicts the same Fermi surface topology and the same 
(tt, tt) nesting [lll |. 

Fig. m shows the calculated band dispersions of a typ- 
ical arsenide material, LaFeAsO (Fig. [2fa)), and FeTe 
(Fig. El^b)). As discussed by Subedi e< al. [ll|, the two 
band structures are similar. The Fe-Sd states are dom- 
inant around the Fermi level while the anion p bands 
depicted by 'fat' bands are located fairly well below the 
Fermi level, and the similar band structure around the 
Fermi level produces the same Fermi surface topology. 
As a result, the same kind of (tt, tt) nesting is obtained 
for FeTe as in the other arsenide materials [ll|. Here we 
focus on the differences found in the electronic structure. 
Firstly, Te-p bands hybridize with Fe-d at around —2.5 
eV along Y-Z line while As-4p is well separated from Fe- 



FIG. 2: (Color online) The calculated band structure for the 
non-magnetic phase of (a) LaFeAsO and (b) FeTe. The As-4p 
and Te-5y) bands are depicted by 'fat' bands. Fermi level is 
set to be zero (horizontal line). 
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FIG. 3: (Color online) Fermi surface (in a-b plane) for FeTe as 
a function of doping level, &: (a) & ~ 0.0, (b) 5 ~ 0.5 electron 
per formula unit. The corner and center of the square unitcell 
corresponds to the F and M point, respectively. 



3d states. Different features along T-Z are also found at 
energies above the Fermi level. In FeTe, there are signif- 
icant band crossings at about -1-0.3 eV. Another notable 
difference exists at around the X point above the Fermi 
level where the parabolic band along X-M-Y is flatten- 
ing at about -1-0.5 eV in FeTe. It is also noted that there 
is no band crossings across the X point in the range of 
+0.5 — 1-1.8 eV. These differences above the Fermi level 
suggest that possibly different Fermi surface topology is 
induced by electron doping as excess Fe atoms appear in 
stoichiometric FeTe host. 

Fig. [3] shows the Fermi surface for Fci+^jTe as a func- 
tion of doping: (a) without doping and (b) doping by 
0.5 electrons per formula unit. Fig.[3]Ja) is in good agree- 
ment with previous calculation |ll| and clearly shows the 
existence of (tt, tt) nesting as in the arsenides. The most 
important feature found in Fig. [Hfb) is the square-type 
topology developed at around F and (tt, 0) points with a 



similar size. It suggests a new nesting at (7r,0), which is 
consistent with the rotated and doubled stripe spin struc- 
ture (Fig. mb)). The doping level hy 6 ^ 0.5 electrons 
would correspond to Fei.oesTe provided all eight valence 
electrons of excess Fe atoms contribute to the change 
in the Fermi level within a simple rigid band approxi- 
mation. According to a recent neutron scattering by Li 
et al, the commensurate doubled stripe (Fig.[ljb)) spin 
structure is realized [l^ in Fei.oesTe which is very close 
to our case. Such agreement assumes that our simplified 
rigid band treatment of doping may indeed capture the 
essential physics of this system. 

The nesting property is further examined by calculat- 
ing the Stoner response function, xo, given by 
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Fig. m^a) shows the imaginary part of xo(qz = 0) for 
the stoichiometric FeTe without doping (no excess Fe) in 
which the (tt, tt) nesting is clearly seen and is consistent 
with its Fermi surface in Fig.[3I^a). The Imxo{<lz = 0) for 
S — 0.5 (equivalent to Fei.oesTe) is shown in Fig. [IJb). 
Note the strong (tt, 0) nesting and the suppressed (tt, tt) 
nesting, which demonstrates the remarkable difference of 
the doped FeTe from the undoped FeTe and Fe arsednies. 
The novel (tt, 0) nesting is derived from the square-type 
Fermi surfaces shown in Fig. [3]Jb) and is driven by the 
electron doping through the excess Fe atoms. It is noted 
that the small amount of excess Fe plays the key role in 
determining the magnetic structure as the result of the 
Fermi surface change. Since this (tt, 0) nesting matches 
with the rotated-doubled spin stripe found in experi- 
ment, our result strongly supports the same spin fluctua- 
tion mechanism for superconductivity in Fe chalcogenides 
as in the Fe arsenides. 

The nesting property of Fei.oesTe at ^ 0, is dif- 
ferent from that at qz = 0. While any notable peak is 
not found in Imx{ciz 7^ 0) for the stoichiometric undoped 
FeTe, the Imxo{(lz 7^ 0) of Fei.oesTe shows rather com- 
plicated features including both (tt, tt) and (tt, 0) peaks 
with reduced intensities. These features are attributed 
to a significant variation of the Fermi surface along the 
.^-direction which is also reflected in the band disper- 
sion in Fig. [21 The doping level of S ^ 0.5 corresponds to 
the Fermi level shifted by about 0.76 eV where it can be 
noted that the FeTe bands along T-Z direction around 
^ H-0.76 eV region are different from those of LaFeAsO 
(Fig.©. 

By further dopings, the Fermi surface topology once 
again changes significantly, and the (tt, 0) nesting dis- 
appears. At the doping level of S ^ 1.1, which would 
correspond in our analysis to Fei.i4iTe with the incom- 
mensurate spin order as observed in the experiment [l9[ , 
the (tt, 0) nesting is largely suppressed. At this doping 
level, square-type structures are no longer found in the 
Fermi surface, and the complicated multiple-peak struc- 
ture is found to be developed in the xo plot. It might be 
responsible for the experimentally observed inconimen- 
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FIG. 4: (Color online) The imaginary part of calculated 
Stoner susceptibility, xo(ciz = 0), for FeTe in arbitrary unit: 
(a) S ~ 0.0, (b) S ~ 0.5 electron per formula unit. 



surate spin orderings [3, [T^. While, at qz = 0, neither 
(tt, tt) nor (tt, 0) peak is found in xo: several prominent 
peaks exist at around F point which are gradually sup- 
pressed along the (tt, tt) line. At q^ = 2tt/3c, notable four 
peaks are found around the (tt, tt) points. This structure 
found at q^ 7^ shows another difference of Fei+a;Te from 
FeAs materials and is originated from the different band 
structures and hybridizations around Z point above the 
Fermi level as shown in Fig. [21 The accounting for the 
excess Fe atoms and the validity of the rigid band ap- 
proximation is a highly non-trivial problem even when 
we tried to simulate the doping by supercell calculations. 
For the present study, however, our excellent agreement 
with experiment demonstrates that our simplified treat- 
ment well describes this material. 

The magnetic properties are summarized in Table [H 
Note that the calculated magnetic moment based on the 
experimental z(Te) is about 2.09 (Xb which is in good 
agreement with the neutron data ~ 1.97/iB. This re- 
sult demonstrates another difference between FeTe and 
Fe arsenides. In the arsenides, it is known that using ex- 
perimental z{As) in the calculation always leads to over- 
estimating the Fe moment j2^, [s^, [HI . For comparison, 
we also present in Table [J the data for parent material, 
LaFeAsO. The numbers in parenthesis are the experi- 
mental ones. The overestimation by LSDA is more than 
a factor of four, which is exceptionally large. The ori- 
gin of this large discrepancy is still under debate Q . If 
the small moment observed in experiment is attributed 
to domain motions as suggested by Mazin and Johannes 
[3^ . the good agreement found in FeTe implies that this 
system is free from such dynamics, which calls the further 
investigation along this line. 

Table U also shows the calculated exchange couplings of 
the Heisenberg spin haniiltonian, H = J'^^ijy S.^ • Sj, 
in FeTe compared to the arsenides. Importantly, with- 
out doping, even on top of the double-stripe spin or- 
dering, our linear response calculations predict the un- 
stable magnetic interactions. While the charge density 
self-consistency is achieved for both single- and double- 
stripe spin order, the spin waves for the undoped FeTe are 
found to be unstable with respect to the spin angle tilt- 
ing. For the single-striped FeTe, Ju becomes AFM and 



4 



System Moment Jia Jib 'ha J26 

double stripe Fei.oesTe 2.09 (1.97") -7.6 -26.5 46.5 -34.9 

(7r,0) FeTe 2.16 -4.2 12.9 -6.2 -15.3 

single stripe FeTe 2m 38.6 21.7 51) ~ 

(tt.tt) LaFeAsO" 1.69 (0.36") 47.4 -6.9 22.4 



is stabilized. In Fei.ogsTe, on the other hand, J2a is the 
strongest while Jia is as small as Ju in the arsenides. It 
is noted that the two FM coupling, Ju and J26, are fairly 
large which might be responsible for the novel doubled 
stripe AFM spin order. Our results suggest the spin wave 
velocities and dispersions for Fe chalcogenides are differ- 
ent from the arsenides, which can be verified by inelastic 
neutron scattering. 

In conclusion, our electronic structure calculations 
show that a small amount of excess Fe atoms existing in 
the Fei+a;Te samples changes the Fermi surface topology 
significantly. As a result, a novel (tt, 0) nesting appears 
at a; « 0.063, and the rotated double-stripe AFM spin 
structure stabilizes. This is different from Fe arsenide 
parent materials. Our finding of the 'missing nesting' ex- 
plains the origin of the spin density wave observed by a 
recent neutron experiment and validates the spin fluctu- 
ation theory of superconductivity for Fe chalcogenides. 
Further doping suppresses the (tt, 0) nesting and pro- 
duces a multi-peak structure in the Stoner susceptibil- 
ity which might be responsible for the incommensurate 
spin order observed in experiments at higher levels of ex- 
cess Fe atoms. The calculated exchange interactions and 
spin moment demonstrate the role of excess Fe atoms in 
stabilizing the magnetic structure, and imply a different 
magnetic behavior of chalcogenide superconductors from 
the arsenides. 
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TABLE I: The calculated Fe moment (in hb) and exchange 
parameters (in meV) for double stripe Fei.oegTe (doped) and 
FeTe (undoped) along with the single stripe FeTe (undoped). 
The results of LaFeAsO are also presented for comparison and 
the experimental moments are given in parenthesis. 

its strength is much larger than J2a which is in a sharp 
disagreement with the single-stripe-stabilized LaFeAsO. 
In the undoped double stripe FeTe, the overall size of ex- 
change interactions is small, and the signs of Jia, Jib, J2a 
do not correspond to the actual ordering. It implies that 
the stoichiometric FeTe {x = 0.0) is hardly stabilized, 
which partly explains the reason that the FeTe sample 
always has some amounts of excess Fe atoms [l^ • 

The spin waves and magnetic interactions are stabi- 
lized in Fei.oesTe being consistent with experiment. In 
arsenides, the exchange interactions are represented by 
two major AFM interactions, Jia and J2, and their 
strengths are in the same range as is seen in the entire 
arsenide family [2^: The Jia ^ 45 meV and J2 ~ 20 
meV while the ferromagnetic Ju is very small. Based 
on these exchange interactions, the striped AFM phase 
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